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For large order, Laguerre polynomials can be approximated by Bessel functions near the origin. This can be 
used to turn many Laguerre identities into corresponding identities for Bessel functions. We will illustrate 
this idea with a number of examples. In particular, we will derive a generalization of a identity due to Sonine, 
which appears to be new. 

Bessel functions are ubiquitous within mathematical physics, and there are a large number of sum and integral 
i^l ", identities known for them. The same is true for Laguerre polynomials. Since a particular limit of the Laguerre 
y—( • polynomials is given by Bessel functions, it is natural to consider whether some of the Laguerre identities could get a 
f^) | new life as Bessel identities. 

What inspired this note, was an identity that showed up in the solution of a random matrix models. Many 
random matrix models are solvable, for finite matrix dimensions N, in terms of orthogonal polynomials^. In general 
the interesting case is when N goes to infinity. In the case of Gaussian ensembles, the solution is in terms of Laguerre 
(or Hermite) polynomials. If one is interested in the spectrum near the origin, one is then naturally lead to the kind 
of calculations we consider here. 

The well known property of Laguerre polynomials we shall make use of, is that, for a real and fixed, 

uniformly for x on some bounded region of C (eg. Theorem 8.1.3 of ref. 0). Using this together with standard asymp- 
_C totic estimates like Stirling's approximation, the passage from Laguerre identities to identities for Bessel functions is 

jrt ' mostly straightforward. This, of course, requires that the occurrences of the Laguerre polynomials are such that one 
can approach asymptotic region. 

A curious phenomenon can occur if Laguerre polynomials with negative are present. When converting sums over 
such polynomials to integrals, it can be necessary to add so-called anomalous terms (as noted in ref. 1). We define 
Laguerre polynomials at negative index by continuation of the usual formula 

£ m(»):=E^rf^ + a V' = £ .,/" )J . M (a + n)(a + n-l)---(a + j + iy. (2) 
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For general a, clearly L^(0) ^ 0. But, from the last expression, we see that the lowest order term of L~"(x), with v 
^_^ • a positive integer, is x v when m > v. This already hints at the fact that something special happens at the negative 
integers. We will also see that this is connected with the behavior of Bessel functions near the origin. 

The idea of using the fundamental relation between Laguerre polynomials and Bessel functions is very general, and 

appears to be relatively unexplored. We seek to illustrate this method, by providing a number of examples, in the 

hope, and with the expectation that, others will find it useful. Some of the identities, like (fTTj) . we prove rigorously, 

while for others, like (|15p . a complete proof is beyond the scope of this paper. We stress that we have checked all the 

5— 1 ■ identities numerically to high precision. 

Let us now consider some concrete examples. In the following, it will be useful to have (Q} in a slightly different 
form. For fixed a£l and x e C, we have 
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lim N-«L? N (^-)=2*r a / 2 x- a J a (Vfx), (3) 
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which holds uniformly for r on a bounded interval of the positive real axis (excluding the origin) . 

I. SONINE LIKE IDENTITIES 

In this section we will derive two identities giving closed expressions for certain integrals of products of two Bessel 
functions. The first one generalizes a formula due to Sonine. 

The starting point is the known sum formula for Laguerre polynomials (eq. 22.12.6 in ref. [5|) 

N 

V L: n (X)L^ N _ m (Y) = L% + f> +1 (X + Y), (4) 



with the abbreviations 

X :=— and Y := Xr=, 
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which we will use from now on. Naturally this holds for arbitrary complex x, y, a, /3, but let us for the moment restrict 
to a,P > — 1. We can now let N —¥ oo on both sides, which will replace the polynomials by Bessel functions. The 
intuition is that, using ([3]) we replace the polynomials with Bessel function, and that the LHS turns into an integral 
over r. 

This intuition turns out to be too naive in some cases, so let us work out the details. Fix x, y =£ 0, and e > 0, and 
set r = m/N. First we choose M such that 

\N- a -n^(X)L ff N _ m (Y) - 2«+^«/ 2 (l - r)^ 2 x- a y-f l J a (^x)MVT^y)\ < e, 
when M < m < N — M. It then follows that 
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Now choose N such that 

N-M 



= M 



< e. 
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To see that the last condition is possible to satisfy because, note that e.g. the first term in the sum grows as 7V^, 
which is dominated by N a+l3+1 , since a + 1 > 0. 
Combining ([S]),© and ([7]) we conclude that 

lim N- a -?- 1 V Lf n {X)L p N _ rn {Y) = 2^x~ a y- fi / drr«/ 2 (l - rf' 2 J a {^x) J (VT^y), 

since e was arbitrary. Taking the corresponding limit of the RHS of (j4]) , we obtain the identity 



a y~P / drr Q / 2 (l - rf' 2 J a (^x) Jp(VT~ry) = 2(x 2 + tf)-lf*+f>+W J a+(i+l (^Ty^). 



(8) 

By complex continuation^, this holds for arbitrary a, ft with Re a, Re /3 > — 1, which ensures convergence of the 
integral. This formula is known as Sonine's second integral 7 , and a proof based on expanding the Bessel functions 
using an integral representation can be found in ref. [8|. Later we will sketch yet another derivation from an Laguerre 
integral identity. 

It is possible to generalize ©. We note that, due to the relation J_ n (x) = (— ) n J n (x), the integral of © is also 
convergent when a and/or j3 are negative integers. However, the identity does not hold as it stands, because we can 
no longer satisfy ([7]). Let us for definiteness consider the case where a is a negative integer. From the discussion after 
([2]), we see that 

lim N- a L^(X) 

N— s-oo 

is only finite when m > —a. We thus have to treat the terms with m < —a separately. Expanding the polynomial in 
X, we find 
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Here A is the forward difference, i.e. A/ m = f m +i — fm, and we have used that (^ l ) = (— ) fc ( ™~t x ). It is 
easy to simplify the difference operator acting on the Laguerre polynomial. Induction on the identity L^ (z) — 
L^(z) - L^_ x (z) yields 

A k L^(z) = L a m k k {z). (10) 

There is no further cancellation among the polynomials. The limit of ([9]) is therefore straightforward, and we find 

°° m=0 j=0 "'■ ^ ^ ' 

Note that these individual terms give contributions of the same order as the rest of the sum. With the difficult terms 
explicitly dealt with, the considerations leading to ((5J) goes through for the remaining terms, and we end up with the 
identity 

\x~ a y-P J drr a / 2 (l-rf/ 2 J a (Vfx)MVT^y) = (x 2 + y 2 )- {a+P+1)/2 J a+ p+i(Vx T T7) 
1 Jo 

-y— "- 1 E 7F (-£) Jo+^+ift/)-*— "- 1 E ?(■&") J ^+i+ 1 ^' (11) 
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which is valid both when Re a > — 1 (with the sum understood as being zero), and when a is a negative integer. The 
same applies to j3. We see that the integral given by the LHS of ([lip is discontinuous at a = — 1 (approaching from 
the Re a > —1 half-plane), with the size of the discontinuity given by the sum term on the RHS. This is related to 
the disappearance of the branch cut of J a (z) along the negative real axis in the z plane at a = — 1. 
Taking a = and /3 = — 1, we obtain 
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dr (1 - r)- 1 ' 2 J (V?x) JiiVT^y) = -JoWx 2 +y 2 ) + J {x), 



which was previously obtained in the context of a random matrix calculation^. This was extended to the series of 
identities with a — (3 + 1 a negative integer in ref. [2|, but we have not found the general result pip in the literature. 
The example just given was special in that we could compute the sum of the terms we had to separate out explicitly 
at finite- N. In general, this will not be possible. To illustrate how to handle this, let us consider the identity (48.21.14 
in ref. |9|) 

where (x) m := x(x+l) ■ ■ ■ (x+m—1) denote the rising factorial (or Pochhammer symbol). For a, j3 > we can, mutatis 
mutandis, repeat the above analysis. The only additional ingredient is that we need to use Stirling's approximation 
for the factorials. The result is 

F(f^i T P"-' ,,2 - 1 <l -rf<W*)MVT=-r l} = rW J ffl + 1) W. (12) 

which appears to be a new identity. 

As before, we can also consider the case where j3 is a negative integer. We again find that the terms with N—m < —f3 
are anomalous, but in this case we will apply a more general technique to calculate the large- N limit. By a rewriting 
similarly to what was done in (0), we have 



V^ {\[a-P\) H V^ 1 (-X)P „ , 
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Lacking the analog of (jlOp , we need a different approach. As shown in the appendix, we can, assuming some regularity 
of the finite- N corrections, calculate the large- N limit. The prescription amounts to taking the limit of the expression 



in the square bracket for general m, and then effectively substitute A m — » N l d r . Carrying this out, we have 

lim N-iM £ { ^~^ )m L a m {-X)L^_ m {X) 
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which lead us to generalize (|12[1 to 
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Convergence of the integral requires Re[a — /?] > and Re/3 > — 1, or that ft is a negative integer. Let us give the 
identity for the first few negative values of ft explicitly 

drr-^ 2 (l-r)- l / 2 I a (V^x)J 1 (VT^x) + I a (x) = ^ [ " + ^ J a (x) , 
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II. ADDITIONAL IDENTITIES 



We briefly mention some other identities that we found interesting. 
Let us first consider the formula (eq. 48.3.13 of ref M) 



m=0 

with v a non-negative integer. The large- N limit leads to 



{-N) m {v + l/2) m r _ 2v (N + v) u _ v 2 
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dr [r(l - r)]- 1/2 J 2 „(2^x) = vr [J.(a;)] 2 
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where we have used that, since v is an integer, J_„ = (— )",/„. If we allow z/ to be a general complex number (with 
Rev > —1/2 to ensure convergence), we cannot deduce from the derivation that the identity should continue to hold. 
Numeric calculations, however, strongly suggest the it remains valid on this enlarged domain. 

The identities we have considered until now have involved sums over the degree of the polynomials, but one can 
also consider sums involving the order. Let us take 



(—) ,n 
° Y J2 ^\-Y m L a N +m {X) = L a N (X + Y) 



m=0 
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as an example. The asymptotic relation © we have used so far is for fixed order, so it is does not bound the error 
of replacing Laguerres with Bessel functions in the sum over m. If we nevertheless proceed with the replacement, we 
are lead to 

^"E^T f J a+m {x) = [x 2 + y 2 ]- a/2 U^ r +y 1 )- (15) 

m=0 ' ^ ' 

We have verified this is a valid identity to high precision numerically. 

Lastly, let us remark that the same idea applies when starting from an integral Laguerre identities. As an example, 
let us take (eq. 22.13.13 of ref. i) 

X dYY a (X - Yy- l L%{Y) = r(/3) r( ^ + "+| ) 1) A-+^L^(X), 

where convergence requires Rea > — 1, Re/3 > —2. Since the integration domain is bounded, taking the large- N limit 
presents no difficulty, and we find 

dyy a+1 (x 2 - yY-'My) = 2^ 1 r(P)x a+ f s J a+P (x). (16) 

It is interesting to note that (8) can also be obtained from an integral identity. Indeed, starting from 10 (2.19.14.1 
in ref. [Ill) 

[ d "°c - #*. (£) 1 (w) - ^ B <° + '" + 1 -" + " + 1 > i »i- +1 {w) • (") 

and letting N, n, m —> oo such that n/N — > x 2 and m/N — > y 2 , we easily recover the previous result. Note, however, 
that (fl~7| only converges for Re a, Re /3 > — 1, so it does not seem possible to obtain the generalization (|llft this way. 

III. CONCLUSION 

We have illustrated a general method for obtaining Bessel identities through a number of examples. The identities 
we have given as examples appear to be new. 
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Appendix: Finite differences in the large- TV limit 

We will give a brief argument that, when the large- N limit is well-behaved, the finite difference operator A effectively 
turns into a differential, in a way that will be made precise. For that purpose, let us consider a family of functions 
Jn{z) : C — > C for each integer N > 0. Assume that we have the asymptotic expansion 






with F, E p analytic. Then we havei^ (I being a fixed integer) 
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For any k > Taylor expansion yields 
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The crucial property of (|A.1|) is that the leading (in A 1 ) coefhcient of Z p is just a derivative of F(z). Using that, for 
any polynomial Pi :— ^2 p N~ p a p l p , we have 



A p P t = p\N- p a p + 0(N- p - L ) 



we get the result 



lim N- p A p f N (z) = |L 
N^>oo or- 1 



F(rz). 
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